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ABSTRACT 


Photoelasticity is an optical method of experimental stress analysis, which yields a 
whole field representation of the principal stress difference and its orientations. With 
advancements in digital image processing techniques applied to photoelasticity it has 
now became possible to get fringe order and isoclinics at every pixel in the domain. 
Usually the retardation calculated is represented in the form of a phase map. Due to 
isoclinic-isochromatic interaction on isochromatics, the phase map has ambiguity in the 
sign of fractional retardation in some zones. The present work discusses different 
methods of improving the phase map at least for particular cases. The effect of slow and 
fast model axis in the calculation of isochromatic parameter is verified experimenfeilly. 

The various difficulties related to the evaluation of isoclinic parameter are 
presented systematically. A new approach to unwrap isoclinic parameter by using 
isochromatic parameter is discussed. Various methods for obtaining the isoclinic 
parameter in plane and circular polariscopes are implemented. A quantitative 
comparative study has been conducted on the performance of these methods on tire 
isoclinic for the problem of ring under diametral compression. The study has showed in 
general the PST methods in plane polariscope are better compared to PST methods 
based on circular polariscope. Among the polarization stepped methods, the method 
proposed by Brown and Sullivan is found to be good. In PST methods in plane 
polariscope, the method proposed by Ramesh and Mangal is found to be good. 
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CHAPTER 1 


INTRODUCTION 


1.1 Introduction 

Photoelasticity deals with the determination of stress fields by use the of the polarized 
light. This method is based on the temporary birefiringence effect. Certain transparent, 
non-crystalline materials like Plexiglass, celluloid, homolite, epoxy etc., are optically 
isotropic when not loaded but become optically anisotropic and display characteristics 
similar to crystals when they are stressed. The birefringence in the material is retained 
only during the application of the loads and disappears when they are removed. 
Measurements of this birefiingence provide data adequate for quantitative determination 
of the corresponding state of stress. The measurements are made with the help of an 
optical instrument called polariscope. A polariscope could be used to find out viz. the 
maximum shear stress i.e., the difference in the principle stresses and the orientation of 
principle stresses. The corresponding contours called isochromatic and isoclinic fiinges. 
The experimentation in photoelasticity involves acquisition of these images and their 
interpretation. The main merit of photoelasticity is its simplicity in experimentation as 
well as in the evaluation of the stress field. 

Photoelasticity directly provides the information of principal stress difference 
and the orientation of the principal stress direction at the point of interest. Using 
these, one can find the normal stress difference and in-plane shear stress by 
invoking equations in mechanics of solids or Mohr’s circle. 
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1.2 Present Work 

In photoelasticity, the phase information is affected by both the difference in 
principal stresses and the orientation of the principle stress direction. This is 
referred to as isochromatic-isoclinic interaction and this affects the evaluation of 
isochromatics and isoclinics differently. Hence in digital photoelasticity the whole 
field evaluation of isoclinic parameter poses certain problems. Use of six-step 
phase shifting technique, which uses the intensity data in spatial domain, in 
general provides the estimation of fringe order with a reasonable accuracy. 
However, the isoclinic parameter is not quite accurate. 

Due to isoclinic-isochromatic interaction on isochromatics, the phase map has 
ambiguity in the sign of fractional retardation in some zones. The present work 
discusses different methods of improving the phase map at least for particular 
cases. A new approach to process image in quadrants is presented here. An 
effective use of taking logical NOT of a pixel value for the improvement of phase 
map is presented. The effect of slow and fast model axis in the calculation of 
isochromatic parameter is verified experimentally. 

The intensity patterns obtained from six-step phase shifting method are 
studied for the problem of disk imder diametral compression. The various 
difficulties related to the evaluation of isoclinic parameter are presented 
systematically. A new approach to unwrap isoclinic parameter by using 
isochromatic parameter is discussed. Various methods for obtaining the isoclinic 
parameter in plane and circular polariscopes are implemented. A quantitative 
comparative study has been conducted on the performance of these methods on the 
isoclinic for the problem of ring under diametral compression. 

1.3 Thesis Layout 

In chapter 2, first six-step phase shifting method is discussed in detail. After that 
new approaches based on processing the images in quadrants to overcome the 
limitations in six-step phase shifting method has been discussed. Chapter 2 also 
discusses in details on, how the orientation of model slow and fast axis influences 
the isoclinic parameter. A new approach of logical NOT to improve the phase map 
is also discussed. 
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Chapter 3 focuses on the extraction of experimental isoclinic data from 
various methods in plane and circular polariscopes. First a study of intensity 
patterns is done in six-step phase shifting technique. The study has been done to 
understand the reasons of numerical inaccuracies involved in the calculation of 
isoclinic parameters. A different way of arriving at isoclinic parameter from the 
isochromatic parameter is studied. 

A quantitative comparative study of various phase shifting methods and 
polarization stepping methods for the calculation of isoclinic parameter in plane 
and circular polariscopes has been done for the problem of ring under diametral 
compression. 

Chapter 4 presents the conclusion and suggestion for the future work. 



CHAPTER 2 


RETARDATION CALCULATION IN DIGITAL 
PHOTOELASTICITY 


2.1 Introduction 

One of the early approaches for retardation calculation over the complete model domain 
was the method of phase shifting. In this, intensity data recorded for a few optical 
arrangements are judiciously used to solve for the isochromatic and isoclinic fringe 
order. Following the lines of phase shifting [1], several vmants of these appeared such 
as polarization stepping [2] and load stepping [3]. In polarization stepping the polarizer 
and analyzer are kept crossed for recording intensity data. Load stepping essentially 
adopts one of the existing phase shifting techniques for the intensity recording. Instead 
of recording intensity data for one load, the data is recorded for small steps of loads on 
either side of the load of interest. The methodologies such as phase shifting, polarization 
stepping or load stepping have their own advantages and inherent limitations, in this 
chapter, the limitations of phase shifting are initially highli^ted. New approaches to 
overcome these limitations, based upon the processing the image in quadrants have been 
proposed. These approaches are validated for the problem of a disk xmder diametral 
compression. 

2.2 Intensity of Light Transmitted for a Generic Arrangement of 
a Circular Polariscope 

Consider the optical arrangement shown in Fig. 2.1 in which the second quarter-wave 
plate and the analyzer are kept at arbitrary positions. The first quarter-wave plate is 
kept at 13 5-deg. 




Fig. 2.1. Generic arrangement of a circular polariscope to employ phase shifting methodology 


Using Jones calculus, the components of light vector along the analyzer axis and 
perpendicular to the analyzer axis (for 135-deg) are obtained as 


E, 1 i 

cos P 

siny0 

l-icaslT] -isin 27 

J 2 

-sin >9 

cos P 

-isin 2;7 l + icos 2;7 


s s s' 

cos i sin— cos 2^ -i sin— sin 2^ 

2 2 . 2 

s s s 

-isin— sm20 cos— + i sin— cos 20 

2 2 2 . 


1 i 
i 1 



imt 


( 2 . 1 ) 


The combination of the quarter-wave plate at 135-deg and the polarizer produces a 
left circularly polarised light. Denoting the intensity of light transmitted as the 
intensity expression is obtained as 

If = ^ + ^ [sin 2(y5 - t]) cos <5 - sin 2(0 - rf) cos - rf) sin 0] (2.2) 

where 4 accounts for the amplitude of light vector and proportionality constant. The 
intensity equation developed show that when the extinction conditions are satisfied, the 
intensity of light transmitted should be zero. However, this never happens in an 
experiment. In order to account for this behaviour, a term h is added to account for the 
stray light/backgroxmd illumination. 

In subsequent discussions, the intensity of light transmitted will be denoted as 
for the i** optical arrangement. The orientation of the optical element that needs to be 



6 


rotated will also appear with a subscript i. Thus, in Eq. (2.2) by adding the term h to 
account for background illumination/stray light reduce to 


h = h + Y + Y[sin2(>9,- - 77 ,- ) cos - sin 2(^ -;/,•) cos 2(y9,- - 77 ,.)sin J] (2.3) 


2.3 Six-step Phase Shifting Aigorithm 

One of the most successful and time-tested approach of phase shifting to evaluate 
isoclinic and isochromatic parameter involves the use of six different optical 
arrangements [1]. The various positions and the respective intensity equations are 
summarized in Table 2.1. The arrangement shown also accounts for the mismatch of 
quarter wave plate error. 

Using the intensity equation listed in Table 2.1, the isoclinic parameter 6 is 
obtained as 





(2.4) 


Table 2.1. Polariscope arrangement and intensity equations for six-step phase shifting method. 


4 

1 

P 

Intensity equation 

37r/4 

7r/4 

7r/2 

A = 4+y(l + cos<5) 

3:7/4 

n/A 

0 

4 = 4+y(1-cos«^) 

3a/4 

0 

0 

=/jj+^(l-sin2^sin<J) 

3;r/4 

72/4 

7r/4 

=/^ +^(l + cos2^sin<^) 

7Z/4 

0 

0 

/j =/^ +— (l + sin2^sin^) 

7r/4 

37^4 


/g =/g +^(l-cos2^sin<5) 


From the equations shown in Table 2.1, the fractional retardation 5 is obtained 
as 



(2.5) 


5^ = tan ^ 


5^ = tan ^ 


(/4-4) 

(/l -/2)cos2^cy 

(/ 5 -/ 3 ) 1 

(/i -I^)sm2e, ) 


= — tan ' 
2 

= — tan~^ 
2 


^ cos 26>(, sin <5^ ^ 
J^cosie^cosS^^ 

I a sin 2^c sin 5^ '' 
J^smie^cosS^ ^ 


( 2 . 6 ) 


In Eqs. (2.4) to (2.6) the subscript c for the Box 5 indicates that the principal values of 
the inverse trigonometric function is referred. In subsequent discussion too the same 
symbolism will be used 

2.4 Advantages and Limitations of the Six-step Phase Shifting 
Method 


The exact calculation of isoclinic angle and fractional relardation is not trivial from Eqs. 
(2,4) -(2.6). Equation (2.4) shows that the isoclinic value is indeterminate at those 
points where 5^ = 0 , it, lit, Zit .... Tlie implication is that the isoclinic contour will not 
be continuous over the domain. When 5^ is not exactly equal to 0 or or 2it etc. but 
very close to these values, 9^ does not become indeterminate but the value of 9^ 
determined will be unreliable and appears as noise in isoclinic fringe when plotted. This 
is known as isochromatic-isoclinic interaction in the calculation of the isoclinic 
parameter. 

Considering that 9^ is calculated correctly by Eq. (2.4), the calculation of 5^ by 
Eqs. (2.5) or (2.6) should result in the same value of 5^ if continuous variables are used 
in the arctangent calculations. In practical applications, this is not the case due to 
intensity digitisation and quantization and only discrete values are used in the arctangent 
calculations. This can lead to error in the calculations. Further, in Eqs. (2.5) & (2.6), 
has low-modulation regions when cos2^(; or sin2^g are small. In these regions, the 
intensity differences — /g) or (/g — 73 ) will also be small. The evaluation of 5^ will 
involve the ratio of two small discrete numbers and will be in considemble error in such 
regions. This is termed as isoclinic-isochromatic interaction in the calculation of the 
fractional retardation. Apart from this, the principal value of 6 *. evaluated by Eq. (2.4) 
lies in the range -^/4 Xo it! \ and an ambiguity exists on whether corresponds to cTj or 
CT, direction over the domain [ 1 ]. 




Fig. 2.3. Whole field representation of fractional retardation, 4 in the form of phase map of a disk under 
diametral compression obtained using Eq.(2.7) with 0c being calculated by various means: a Oc by Eq. 
(2.4) and 4 by Eq. (2.7) b 6c by theory corresponding to the cTi direction and 4 by Eq. (2.7) c phase map 
Fig. 2.3a improved by interactive approach d 6c by Eq. (2.4) then unwrapped by Eqs. (2.8) & (2.9) 
according to quadrant and 4 by Eq. (2.7) 

It is instructive to note that regions of low values of cos26>, coincide with high values of 
sin2^(. and vice versa. Therefore for every point in the region of interest one has two 


estimates of 5^ with different degrees of error depending on the region. Thus, the results 
must be combined appropriately to get the best possible values of For the whole 
domain it is obtained as [4] 

^ (/g -/3)sin2^e +(A -/6)cos2^c ) ^ 


5^ = tan ' 


(2.7) 


Unlike Eqs. (2.5) or (2.6), Eq. (2.7) always guarantees a high modulation over 


the field. 




Figure 2.2 shows the sequence of six images recorded for a disk under diametral 
compression (dia 60.10 mm, thickness 6.28 mm and load 773.0 N). The material stress 
fringe value is calculated by a digital technique and is found to be 11.9973 N/mm/fr. 
Figure 2.3a shows the phase map obtained by Eq. (2.7). on scrutiny one can see that 
near the load application points the phase map is not good. This is not just due to high 
stress gradient but is due to the interaction of isoclinic on the isochromatic evaluation. 
The principal values of 0^ of Eq. (2.4) lies in the region -^/4 ^ 0c ^ ^4 whereas, 
physically, principal stress direction lies in the range -^/2 <0< ^2 and this introduces 
an ambiguity on 5 evaluation. Figure 2.3b shows the phase map obtained from 
experimental images using the theoretically calculated 0 indicating the ctj direction 
throughout the domain. The phase map is uniformly good over the domain. 

Though, phase shifting can provide 5 over the whole domain, certain level of 
human interaction is needed to correctly identify its sign for unwrapping. This is one of 
the serious limitation of the technique and could be resolved only when isoclinic is 
unwrapped correctly over the domain corresponding to either ctj or <^2 direction 
uniformly over the domain [1,5]. This problem is a current area of research and 
automatic solution is still elusive. Figure 2.3c shows the phase map of 2.3a improved by 
the interactive approach developed by Ramesh and Mangal [1]. The approach is a 
general one and is applicable for any problem and one requires the information of fringe 
gradient in the bad zone for the phase map improvement. 

The basic advantage of phase shifting is that inverse tangent functions are used 
for ^calculation. The unwrapping of (^becomes a lot more simpler and straight forward. 

2.5 New Approaches for Plotting the Phase Map 

The phase maps are influenced by the value of isoclinic parameter, particularly by it’s 
sign. The idea is to study the effect of various methods of evaluation of the isoclioic 
parameter value on the resultant phase map. Different ways of obtaining the isoclinics 
such as, based upon the circularity of trigonometric functions, use of absolute value of 
intensity difference, use of only positive values of twotheta are discussed here. The 
purpose is to explore various possibilities in plotting the phase map, which may work 
satisfactorily at least for particular cases. Even when the method works on the 
theoretically generated images, it may have problems while processing experimentally 
recorded images. 




Fig. 2.4. Sequence of six images generated theoretically for a disk under diametral compression (dia 
60.10 mm, thickness 6.28 mm and load 773.0 N) for the optical position listed in Table 2.1 






The example problem chosen are disk imder diametral compression (both theoretical 
and experimental) and ring under diametral compression (theoretical). Figure 2.4 shows 
the theoretically generated image for disk under diametral compression and Figure 2.5 
shows for the ring mder diametral compression 

2.5.1 Quadrant Approach 

In this the phase map is plotted by dividing the image into four quadrants All these 
approaches are based upon Eq. (2.7). 

2.5.1. 1 Using the circularity of tan function 

In this method the image is divided into four quadrants as shown in Fig. 2.6. 



Fig. 2.6. Relative positions of four quadrants 


In the first and third quadrant after calculating the theta, the following conditions are 
applied. 


0 = 


e^+Tvli 


if (^ > 0) 

if(^ < 0) 


In the second and fourth quadrant the condition applied is 


if(0<o) 

if(^>0) 


( 2 . 8 ) 


(2.9) 


Changing ftie value of theta by nil will change the values of sin20, cos2^from positive 
to negative. Figure 2.7a shows the phase map obtained by employing the above 
approach for the disk under diametral compression recorded experimentally as in Fig. 
2.2. Figure 2.7b shows the phase map calculated based on theoretically generated 
images as in Fig. 2.4. The phase map is good for the theoretically generated images. For 
the experimentally recorded images the phase map is not good at points very close to 



the vertical and horizontal diameter. This is because of non-symmetric intensity 
patterns. Also the loading direction is not perfectly vertical in the experiment. The 
fractional fringe order for the experimentally recorded images of disk are shown in the 
Fig. 2.8. For the problem of ring under diametral compression, the phase map obtained 
by using theoretically generated images is shown in Fig. 2.7c. This shows that though 
the method works for the disk under diametral compression, it fails for the problem of 
ring under diametral compression. 



Fig 2.7 Phase maps generated by using circularity of tan function a for experimentally recorded images of 
disk under diametral compression b for theoretically generated images of disk under diametral 
compression c for theoretically generated images of ring under diametral compression 
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Fig. 2.8. comparison of fractional fringe order variation for 
(method 1) a along diametral line b along rl2 line c along 3, 




2.5.1 .2 Using the absolute values of intensity difference. 


In this method, the absolute values of - /g) and (/g - I^) are used in Eq. (2.4) to 
calculate 9 . Since both the values involved are positive, the resulting value of ^ will be 
positive only. This positive value of theta is plotted in the first quadrant. In the second 
quadrant same value with negative sign (i.e. with -0 ) is plotted. The results for the disk 
under diametral compression for the experimentally recorded images is shown in Fig. 
2.9a and for the theoreticsilly generated images in Fig 2.9b . In this case, there are some 
noise points in the theoretically calculated phase map. Also areas which are very close 
to the loading point are bad for both experimental and theoretical images. The fractional 
fringe order for the theoretical and experimental images of disk are shown in Fig. 2.10. 
For the ring under diametral compression, the theoretical image is shown in Fig. 2.9c 
This shows that though the method works for the disk under diametral compression, it 
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Fig. 2.10. Comparison of fractional fringe order variation for the mediod of using absolute value of 
intensity difference (method 2) a along diamefral line b along rJ2 line c along 3r/4 line 
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2.5.2 Using Only Positive Value of Twotheta 


In this approach only positive value of 2^ is used in Eq. (2.7). The condition 
applied is, if 2^ is less than zero then ;ris added. The result is shown in Fig. 2.11a for 
the experimentally recorded images and in Fig 2.11b for the theoretically recorded 
images. The phase map is good for the disk in first and third quadrant while it is bad in 
second and fourth quadrant for theoretically and experimentally generated images. The 
phase map calculated based on experimentally recorded images is bad in the regions 
close to vertical and horizontal diameter. This happens because of non-symmetric 
intensity patterns. Also the loading direction is not perfectly vertical in the experiment. 
The firactional fiinge order plot for the experimentally recorded disk is shown in Fig. 
2.12. Fig 2.11c shows the result for theoretically generated images of ring under 
diametral compression. This shows that technique applied to the ring does not give the 
expected results. 



Fig 2.11 Phase maps generated by using positive value of twotheta a for experimentally recorded images 
of disk under diametral compression b for theoretically generated images of disk under diametral 
compression c for theoretically generated images of ring under diametral compression 












2.6 Effect of Swapping Images on Phase Map 

This study is mainly based on the Eq. (2.7). This study is conducted on theoretical and 
experimental images of disk under diametral compression ( load 773N, dia 60. 1 0 mm 
For =11.9973 N/mm/fir). The equations given for six different arrangements in six-step 
phase shifting algorithm are given in table 2.1. These equations are derived by 
considering 6 for the slow axis. Same exercise can be repeated by considering 0for fast 
axis. The results are given in Table 2.2. It can be seen that I\ and 1% are insensitive to 
slow or fast axis but for remaining equations, the signs of sin2fein^ and cos2fein(? are 
changed. 


Table 2.2. Polariscope arrangement and intensity equations for six-step phase shifting method with 
reference to 0at fast axis 
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Let us first study the effect of the loading direction on the intensity patterns. The 
effect of loading direction can be explained by a schematic shown in Figs 2.13a & b. 
Figure 2.13a shows a disc is vertically loaded while Fig 2.13b shows disc loaded on 
horizontal diameter. The orientation of slow (S) and fast (F) axis is shown in both the 
cases. It can be seen that as the loading direction is changed (model along with loading 
is rotated ccw by 90°) the orientation of the slow and fast axis is also changed spatially. 
When loading direction is changed by 90° the fast and slow axis are also rotated by 90° 
spatially. Note that the relative direction with respect to the loading direction is 
constant. But with respect to the light source the positions of both the axis are swapped. 
The tick marked portions shown in two figures will have similar gradient in the phase 
map. A stress firozen disk is used to record images with the loading axis kept vertical 
and loading axis kept horizontal. Fig 2.14a Sc b show the set of images of the stress 
frozen disk kept vertically and horizontally .The Fig 2.15a & b shows the experimental 
phase map obtained with vertical and horizontal loading directions . Referring to Table 
2.2 and Fig. 2.11 one can see that change of loading direction can be equated to the 
swapping the images (I 4 , /s) and (J^ , Ij)- The phase map obtained by swapping the 

images in vertical loading is shown in Fig. 2.15c 



Fig 2.14 a Sequence of six images recorded experimentally for a stress-frozen disk with loading direction 
kept vertically 



Fig 2.14 b Sequence of six images recorded experimentally for a stress-frozen disk with loading direction 
kept horizonally 



Fig 2.15 Phase maps obtained for the experimentally recorded images of stress-frozen disk a For vertical 
loading direction b for horizontal loading direction, c swapping the images ( h , h) &{ for vertical 
loading 


2.7 Obtaining Phase Map by Using iogical NOT 

The Phase map obtained by Eq. (2.7) for theoretically generated images is shown in Fig 
2.16a. The images are processed in the sequence given in Table 2.1. This phase map is 
swapped beyond 45° isoclinic. Fig 2.16b shows a phase map obtained by swappmg the 
images (A , A) and (A , A). This phase map is correct beyond 45° isoclinic. Theoretically 
it is possible to patch the regions beyond 45 isoclinic to obtain a correct phase map. 
This can be done by correctly identifying the areas beyond 45 isoclimc then reading it 
from Fig 2.16b and putting it in the respective position in the Fig. 2.16a. The other way 
of doing it is that after identifying the swapped areas in Fig 2.16a replace the pixel by 
it’s logical NOT. The result is shown in Fig 2.16c. This is done by using a standard 


software Matrox Inspector”. It has a provision to select an arbitrary region of interest 
in the given image and process that part only. The result for the experimental images is 
shown in Fig 2.17a & b. ( load 773N, dia 60.10 mm, F<, =1 1.9973 N/mm/fr). Thus it is 
possible to correct the phase maps by using the logical NOT. 



Fig 2.16 Ph^e maps obtained for theoretically generated images of disk imder diametral compression, 
a without swapping the images b by swapping the images ( I 3 , 1 5 ) <£( /^ , /«) c phase map corrected by 
taking the NOT of pixel value in Fig 2.16 a 



Fig 2.17 Phase maps for experimentally recorded images of disk shown in Fig 2.2 a phase map as per the 
image sequence given in Table 2.1 b improved phase map by taking NOT of pixel value 


2.8 Closure 

Various approaches to obtain a correct phase map in the domain are discussed with their 
limitations in this chapter. The role of arbitrary orientation of model fast or slow axis is 
discussed here. The effect of swapping the sequence of intensity maps has been 
discussed. The use of logical NOT for correcting the phase maps is proposed. 

The main reason for not having a correct phase map all over the domain is that 
the isoclinic parameter obtained is unwrapped. If it is possible to unwrap the isoclinic 
parameter in the range -7iJ2 to +Td 2 , then a correct phase map all over the domain can be 
obtained. Next chapter deals with the comparison of isoclinic parameter in various PST 
algorithms. Also an approach to obtain isoclinic parameter by using the isochromatic 
parameter is discussed. 



CHAPTER 3 


ISOCLINIC PARAMETER CALCULATION IN DIGITAL 

PHOTOELASTICITY 


3.1 Introduction 

The evaluation of isoclinic parameter is comparatively easier while determining it 
manually but the automated evaluation has several problems. The problems concerning 
the whole field isoclinic parameter evaluation can be roughly divided into the following 
parts, 

1. Numerical inaccuracies. 

2. Ambiguity in experimentally evaluating the isoclinic parameter. 

The problem of getting numerically correct value depends on the accurate recording and 
processing of the intensity information. In evaluating this generally arctangent fimctions 
are used in calculations. In practical applications due to intensity digitization and 
quantization, only discrete values are used in the calculations. Also in the regions close 
to the fringe skeleton, the error increases because of the ratios of small discrete values. 
Also there is a significant difference in the intensity maps of the images in theoretical 
and experimental conditions. These reasons lead to an error in the calculations. Besides 
this, in most of the equations the values returned by the arctangent functions are in the 
range oi-TdA to +jdA only, which creates the need for unwrapping. It is to be noted that 
calculations involve inverse trigonometric operations and thus one has a multi-valued 
solutions, this poses a problem of correctly evaluating the isoclinic parameter. Besides 
the numerical inaccuracies, there is an ambiguity involved in evaluating the isoclinic 
parameter experimentally. As it has aheady been discussed in die previous chapter that 
model slow and fast axis influence the intensity patterns. In a model, a priori one does 
not know whether the fast or slow axis lies the range -jdA to +7dA. This ambiguity leads 
to unsuccessfiil determination of 4- Thus the goal of isoclinic parameter evaluation 



shoiild be such that it either gives <j\ or 02 direction uniformly all over the domain. This 
problem could be solved if theta is unwrapped in the range of —nl 2 to +7d2 

In this chapter we shall discuss the reasons of numerical inaccuracies based upon the 
study of intensity pattern. A comparative study has also been carried out on the 
evaluation of isoclinics by different algorithms. 

3.2 Study of Intensity Pattern Over the Domain 

The intensity information represents the properties of the stress field and the orientation 
of fast or slow axis of the model plays a role on the intensity values. The purpose of this 
study is to explore the possibihty of identifying the orientation of fast and slow axis of 
the model on purely the intensity values. For the six step phase shifting algorithm, the 
intensity information considered includes intensities from 4 , 4 & /j to h, the intensity 
differences (4 - 4) and (4 - 4). Following set of equations is derived firom the 
equations given in Table 2.1. 


/i-/2=4cos^ (3.1) 

4 - 4 = 4 sin 29 sin 5 (3.2) 

4-4 =4 cos 2^ sin J (3.3) 

Hh -hf -hf ' (3-“) 

h=(.h + h-h)l'^ ( 3 . 5 ) 


3.2.1 Study of Intensity Variation in Theoretical Images 

This study is done on theoretically generated images shown in Fig. 2.4. The plots 
obtained for the diametral line (y/R = 0.0) are shown in Figs 3.1, 3.2 and 3.3. Figure 3.1 
shows the intensity plots for images I\ to 4 for the diametral line. Figure 3.2 a, b & c 
show the variations of intensity differences (4 — 4 ), (4 ~ 4 ) 3^^ (4 — 4 )- Finally Fig. 
3.3 represents the plots of calculated 4 ? 4 ^nd sui2ftin^, cos2ftiin5 Sumlar results are 
obtained for a line close to the stress concentration zone (y/R = 0.75). They are shown 
in Figs. 3.4, 3.5 and 3.6 respectively. 
















Intensity difference 



Fig. 3.5. Theoretical plot for line at y/R = 0.75 a ( Ii- 12 ) b for ( Ij- 13 ) c for ( Ir- le) 








Fig. 3.6. Theoretical intensity plot for line at y/R = 0.75 a !», Ib b for ( I 5 - 13 )/ 1, c for ( I 4 - !« )/ 

It can be observed in Fig‘3.1 and Fig 3.4 that the intensity maps are symmetric about the 
loading axis. Also the variation is smooth and periodic. As we move towards the 
loading point, the frequency is increased. The amplitude is constant in case of h and h. 
In case of h and k the amplitudes are decreasing towards the loading line and in the 
case of 4 and h the amplitudes are increasing towards the loading line. Same trend is 
reflected in case of intensity differences (/i - 4), {h - h) and (4 - 4 ). But the only 
things which remain unchanged are 4 and 4. Theoretically the value of 4 and /b remain 
constant at 255 and 0 respectively all over the domain. 
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3.2.2 Study of Intensity Pattern In Experimental Images 

This study is done on experimentally recorded images shown in Fig. 2.2. The plots 
obtained for the diametral line (y/R = 0.0) are shown in Fig. 3.7, 3.8 and 3.9. Figure 3.7 
shows the intensity plots for images I\ to Is for the diametral line. Figure 3.8 a, b & c 
show the variations of intensity differences (/i - li), (1$ - h) and (4 - 4). Finally Fig. 
3.9 represents plots of calculated values of 4 4 and sin2ftin(5, cos26fein«5: Similar 
results are obtained for a line close to the stress concentration zone (y/R = 0.75). They 
are shown in Fig. 3.10, 3.11 and 3.12 respectively. It can be observed in Fig 3.10 and 
Fig 3.1 1 that the intensity maps are non-symmetric about the loading axis. Even though 
the variation is smooth, it is not exactly periodic. As we move towards the loading 
point, the frequency is increased. The amplitude also is not constant. In case of h and 4 
the amplitudes are decreasing towards the loading line and in the case of h and 4 the 
amplitudes are increasing towards the loading line. Same trend is reflected in case of 
intensity differences (4 -^ It), (I$ - h) and (4 - 4). It can be seen in Figs. 3.9a & 3.12a 
that 4 and 4 are not constant on the given line of interest. Here 4 and 4* do not remain 
constant all over the domain. Since 4 is a ftmction of time, 4 also becomes a function of 
time. Both the parameters change with the position also. Because of this it is not 
possible to remove the effect of 4 and h completely. This is one of the reasons affecting 
the accuracy of the isoclinic parameter. 





















3.3 Evaluation of Isoclinic Parameter by using Retardation 
Parameter 

Generally isoclinic parameter is calculated first and then it is used to calculate the 
isochromatic parameter. Correct determination of isochromatic parameter depends on 
whether the isoclinic value determined represents uniformly either oi or 02 over the 
domain. In the previous section an attempt was made to identify the cri or 02 direction 
based on the intensity values, typically I 3 . Referring to Eq. (3.2) & (3.3) the sign of 
R.H.S. depends on both 0 and d Hence no general conclusion was possible. On the 
other hand if d is known correctly, it is possible to evaluate d Recently Dhiraj [6] has 
proposed a method to calculate d unambiguously over the domain. In this section the 
possibility of arriving at the correctly imwrapped value of 0 by using d is studied 
theoretically. The results are obtained on theoretically generated images. The 
expression derived is based on die six step phase shifting algorithm. 

From Eq. (3.2) 


sin2^= 

I^smd 

(3.6) 

2sin^cos^= ^ ^ 

I^smo 

(3.7) 

cos 20 = 

sin^ 

(3.8) 

2cos^^-l= / . ^ 

IgSmd 

(3.9) 

2cos^«= +1 

I^smd 

(3-10) 

/„sind 

(3.11) 


Dividing Eq (3.7) by Eq. (3.1 1) 


tan (9 = 


(A-is) 

(A - -^ 6 ) + 4 sin ^ 


(3.12) 


This equation gives the value of theta in the range —tzI2 to +7c/2. Fig 3.13 shows the 
theoretically generated images for the disk. Figure 3.14 shows the phase map obtained 
by quadrant approach.. The isoclinic plot obtained for theta in case of theoretically 
generated disk by using Eq 2.4 are shown in Fig. 3.15a- While Fig 3.15b. shows 
isoclinic plots obtained by using Eq. 3.12. Fig 3.16 shows the results obtained on line 
y/R = 0.9375. It can be seen that the theta obtained by new formula puts the mwrapped 
value of theta. It very closely follows the line for theoretical theta over the line of 
interest. But it also contains spikes where the results by Eq. 2.4 shows the spikes. In this 

case the spikes are higher. This will happen where 8 = Q, tc, 2m . Thus 

theoretically it is possible to imwrapp theta by using 5 . Theta obtained so is continuous. 
But experimentally this may not be the case, since the value of 5 obtained all over the 
domain by load stepping may have certain noise points or some errors. These things 
may affect the theta calculation. Secondly experimentally obtained theta may have more 
discontinuities in the domain. Hence one needs to check for such problems. This 
demands a post-processing of obtained results from the user. 




Fig 3.14 Phase map obtained by quadrant approach for the theoretically generated images given in 
Fig. 3.13 


Fig 3.15 Isoclinic plots for the theoretically generated images of disk under diametral compression a by 

method of Patterson & Wang ( Eq. 3.13) b by new approach ( Eq. 3.12 ) 


F&tterson&yAfeng 

■ By tan ( theta )©<pression 



Fig. 3.16 Plot of theoretical theta by Patterson & Wang, theoretical theta by modified expression, 
theta by theory. 

ZA Evaluation of Isoclinic Parameter in Circular Polariscope 

For the study of evaluation of isoclinic parameter in circular polariscope, six step phase 
shifting algorithms proposed by Patterson & Wang and Ajovalasit et. al. are considered. 

The study is done for the problem of ring imder diametral compression and for the 
problem of disk under diametral compression. The study is done on two different loads. 

3.4.1 Algorithm by Patterson and Wang. 

Arrangements suggested by Patterson and Wang [1] are summarised in Table 3.1. 


Table 3.1, Optical arrangements for Patterson and Wang. 
Patterson and Wang 
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Using the equations listed in Table 3.1, the isoclinic parameter is obtained as 

forsin^T^O (3.13) 


9 =— tan ^ 
2 



= itan-> 

/aSin0sin20^ 

V4--^6> 

2 

sin 0 cos 20^ 


Use of atan() function is recommended to calculate theta from the above equation. The 
experimentally recorded images are shown in Fig. 3.17 


3.4.2 Algorithm by Ajovaiasit et. al. 

Ajovalasit et al [1]. proposed the use of both left and right circularly polarised 
light to arrive at an optical arrangement that minimises the error due to mismatch 
of quarter-wave plate error. The optical arrangements proposed by them and the 
intensity equations are shown in Table 3.2. Arrangements suggested by Ajovalasit et 
al. are summarised in Table 3.2. 


Table 3.2. Optical arrangements for Ajovalasit et. al. 


Ajovalasit et. al. 
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Using the equations listed in Table 3.2, the isoclinic parameter is obtained as 



/„sia^sin2g 
sin ^ cos 20 


for sin0 ^ 0 


(3.14) 


use of atan() function is recommended for the solution of above equation.The 

experimentally recorded images are shown in Fig. 3.18 




Fig 3.17 Sequence of six images recorded experimentally for the method by Patterson & Wang, “ ^ 
the sequence given in Table 3.1 for the ring under diametral compression ( load o03 N, dia. 80 mm, 

thickness 5.05 mm, = 1 1.2329 N/mm/fringe ) 



Fig 3.18 Sequence of six images recorded experimentally for the method by Ajovalasit et al. as per the 
sequence given in Table 3.2 for the ring under diametral compression ( load 503 N, dia. 80 mm, thickness 
5.05 mm, F,, = 1 1 .2329 N/mm/fringe ) 


Fig 3.19 Isoclinic plots for the experimentally recorded images of ring under diametral compression a by 
the method of Patterson & Wang b by the method of Ajovalasit et al. 



Fig 3.20 Isoclinic plots for the experimentally recorded images of disk under diametral compression for a 
low load (load 84 N, Dia. 60 mm, thickness 5.05 mm, F<j= 11.2329 N/mm/fringe ) a by the method of 
Patterson & Wang b by the method of Ajovalasit et al. 

Figures 3.19 a & b show the isoclinic plot for ring from Patterson & Wang and 
Ajovalasit et. al. respectively. The plots show that the isoclinics are discontinuous. 
In some areas they are very poor. Such data is not useful for further analysis. 
Figure 3.20 shows the isoclinic plots obtained for a problem of disk at low load. 
Here the isoclinics are good as compared to the isoclinics obtained at higher loads. 
Thus comparatively better results can be obtained by using a lower load. But this 
may not possible for every problem, particularly for a slice cut from a stress frozen 
disk. Hence, study was carried out in plane polariscope to evaluate the isoclinics at 

higher loads. 




3.5 Intensity of Light Transmitted for a Generic Arrangement of a Plane 
Polariscope 


y 



Fig. 3.21. Generic arrangement of a plane polariscope to employ phase shifting methodology 


Figure 3.21 shows a photoelastic specimen kept in a plane polariscope with the 
polarizer and the analyser kept at arbitrary angles of a. and yd respectively. Using 
Jones calculus, the components of light vector along the analyser axis and 
perpendicular to the analyser axis are obtained as 



cos P sin p 
- sin P cos P 


cos a 


S .. 6 _ 

cos — ism— cos 2d 
2 2 

-i sin— sin 2d 
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5 1 

-isin— sin2d 
2 

5 .. 5 . _ 

cos — + ism— cos 2d 
2 2 


(3.15) 


the intensity of light transmitted /p (the subscript ‘p’ denotes that the incident light 
is plane polarised) is obtained as 


/p =/Jcos''-cos2(yd-a) + sin^|cos2(yd + a-2d) (3.16) 

where accounts for the amplitude of incident light vector and the proportionality 
constant. It is interesting to note that Eq. (3.14) is not altered when the model slow 

axis orientation is changed from d to d+ ;r/2. 
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3.6 Phase Shifting Techniques Based on Plane Polariscope 

Similar to the PST in circular polariscope, different algorithms are developed in plane 
polariscope also. Here a comparative study of some of the main algorithms in plane 
polariscopes is presented. The algorithms considered in plane polariscope are. 

1 Ramesh and Mangal 

2 Brown and Sullivan 

3 Chen and Lin 

4 Dupre etal. 

5 Sarma et al. 

6 Kihara 

In the plane polariscope methods considered, the methods of Brown & sullivan and 
Chen & Lin are actually called as polarization stepping methods while the remaining 
methods are called as phase shifting methods. T his study is done under similar 
conditions to obftdn a comparative statement. This comparison is done for the 
evaluation of isoclinic parameter only. Initially each algorithm is introduced with 
respect to the opticd arrangement and the equation for isoclinic parameter. 
Experimentally recorded images are shown for each optical arrangement m all the 
algorithms. Finally a comparative study of results obtained is done. 


3.6.1 Algorithm by Ramesh and Mangal 

In this, the polarizer and analyzer are kept at arbitrary positions. The intensity of 
light transmitted is given by Eq. (3.16). If the term /j, is added to account for stray 
light/background illumination, one gets 




cos^— cos^(y9,- -a/) + sin^— cos^(^i +a,- - 26 ) 
2 ^ ■ 


(3.17) 


If one takes all possible combinations of a and ^ in steps of then only six 
intensity equations turn out to be unique. The arrangement suggested by Ramesh 
and Mangal [1] is given in Table 3.3. 



Table 3.3. Table showing the intensity equations for different 
orientations of the polarizer and analyzer. 


a 

p 

Intensity equation 



0 

nil 

/j = /b + /a sin^— sin^20 


^/4 

nil 


1 • 2^ • Ar% 

1-sin ~sm40 
2 


0 

nl4 

^3=4+y 

iS 

1 + sin — sm4^ 
2 


/r/4 

3nl4 


Sin —cos 2^ 

2 


tvIA 

nl4 


cos^— + sin^— sin^20 

2 2 

0 

0 


■yS . O S 0 

COS — + sm — cos 20 
_ 2 2 


From Table 3.3, using only the first four intensity equations, the isoclinic 
parameter is calculated as, 


^ =itan~^ 
® 4 


\U-hj 


= itan-> 
4 


i5 

/.sin —sin 4^ 
^ 2 

L sin —cos 40 
2 


.2^ 


forsin*^— 5*0 (3.18) 


Use of atan2 function is recommended for evaluating 6^ from Eq. (3.18). Figure 3.22 
shows the experimentally recorded images for the ring. 




Fig 3.22 Sequence of six images recorded experimentally for the method by Ramesh & Mangal, 
as per the sequence given in Table 3.3 for the ring under diametral compression ( load 503 N, 

Dia. 80 mm, thickness 5.05 mm, Ftr = 1 1 .2329 N/mm/fringe ) 



3.6.2 Algorithm by Brown and Sullivan 


In this polarizer and analyzer are kept crossed and four images are obtained at step of 
7r/8 radians. The polarisation steps used by Brown and Sullivan [1] and the 
appropriate intensity equations are summarised in Table 3.4. 


Table 3.4. Polariscope arrangements and intensity equations 
for polarization stepping schemes. 


a 

Intensity equation 

Brown and Sullivan 
(Polarizer & Analyzer crossed; 

0 

/i = /jj + sin^~sin^2^ 

m 

h =4 + — sin ^^(1- sin 40) 

^r/4 

4 = 4 + 4 sin^^cos^20 

Inti 

4 =4 + — sin^— (l + sin40) 


The isoclinic parameter is obtained as 


9^ = — tan“* 
" 4 


r 


\h 




= — tan * 
4 


/- sin^— sin40 
* 2 

/, sin^|^cos46> 


. 2 ^ 


for sin — *■ 0 
2 


(3.19) 




use of atan2() function is recommended for evaluating above equation. Figure 3.23 
shows the xperimentally recorded images for the ring. 


Fig 3.23 Sequence of four images recorded experimentally for die method by Brown & Sullivan, as 
per the sequence given in Table 3.4 for the ring under diametral compression ( load 503 N, dia. 80 mm, 
thickness 5.05 mm, F„ = 1 1.2329 N/mm/fringe ) 



3.6.3 Algorithm by Chen and Lin 


Chen and Lin [1] proposed another polarisation-stepped scheme. In this initially 
three bright field images in plane polariscope arrangement were obtained. From 
these bright field images, dark-field images are constructed using trigonometric 
identity and the concept of image division. For performing image division, the 
intensity corresponding to unloaded images are also to be recorded. Appropriate 
intensity equations are summarised in Table 3.5. 


Table 3.5. Polariscope arrangements and intensity equations 
for polarization stepping schemes. 


Chen & Lin (Polarizer & Analyzer parallel; ) 


0 




l-sin^2^sin^— 

2 


^/8 J - T 

1- — (sm2^ - COS 20)^ sin^— 

^/4 T — T 

-‘3 ~-‘a 

l-cos^20sin^— 

L 2j 



The intensity equation corresponding to dark-field is obtained by 


(3.20) 

where, is the intensity recorded using an unloaded model. The isoclinic 
parameter is then obtained as 




+ ■^3 ~ 2/2 

h~h y 


(3.21) 


Use of atan2() function is recommended to evaluate 6 ^ from the above equation. 
Figure 3.24 shows the experimentally recorded images for the ring for this arrangement. 




Fig 3.24 Six images recorded experimentally for the method by Chen & Lin, ^ per Table 3.5 
ring under diametral compression ( load 503 N, dia. 80 mm, thickness 5.05 mm, F, - 


for the 
11.2329 


N/mm/fringe ) 



The construction of normalized dark field images from the bright field images 
is an important step in this method. Theoretically is constant at 255 all over the 
domain, and this is always greater than or equal to the loaded bright field intensity 
in the domain. Hence theoretically, generation of dark field images is good. 
Experimentally this is not the case, at some places in the domain the unloaded 
model intensity becomes less than the loaded model. This happens mainly because 
of varying intensity of light source with time. So in the image division the ratio of 
loaded intensity to unloaded one becomes greater than one. Thus the value of I\ in 
Eq. 3.20 becomes — ve at those points in the domain and it is not possible to get the 
correct normalized dark field images. This is rectified by introducing an 
adjustment factor k, so Eq. 3.18 is modified as 






(3.22) 


typically the value of k for the experimental set up used lies in the range 0.5 to 0.7 but it 
may vary in different experimental set up. Theoretically generated normalized dark field 
images are shown in Fig. 3.25. Figure 3.26 shows the experimentally generated dark 
field images without any adjustment (k = 1) and Fig. 3.27 shows the experimentally 
generated dark field images with adjustment (k = 0.65) for the experimentally recorded 
images of ring under diametral compression. 




Fig 3.25 Normalized dark field images generated by Eq. (3,20 ) using theoretically generated images as 
per the sequence given in Table 3.5. a for 0 deg. b for 22.5 deg c for 45 deg 




Fig 3.26 Normalized dark field images generated without adjustment ( Eq. 3.20 ) using experimentally 
recorded images as per the sequence given in Table 3.5. a for 0 deg. b for 22.5 deg c for 45 deg 


Fig 3,27 Normalized dark field images generated with adjustment ( Eq. 3.22, k = 0.60 ) 
using experimentally recorded images as per the sequence given in Table 3.5. a for 0 deg. b for 22.5 deg 
c for 45 deg 


3.6.4 Algorithm by Dupre et al. 



The isoclinic parameter is calculated as 


^ 1 -if 2/3 -/i] 1 


sin —sin 20 ^ 

-1 2 forsin^-?tO 

sin^— cos2^ ^ 

I 2 j 


(3.23) 


use of atan2() function is recommended to evaluate this equation. Figure 3.28 shows a 
set of experimentally recorded images. 




Fig 3.28 Sequence of Four images recorded experimentally for die method by Dupre et al, as given in 
Table 3.6 for the ring under diametral compression ( load 503 N, dia. 80 mm, diidkness 5.05 mm, F^ = 
1 1 .2329 N/mm/fringe ) 


3.6.5 Algorithm by Sarma et al. 

Arrangements suggested by Sarma et al. [1] are summarised in Table 3.7. 


Table 3.7* Optical arrangements for Sarma et al. 






O!) 


The expression for the determination of isoclinic parameter by Sarma et al. method 


IS 


6 . = — tan“* 


27, 


vA +^2 -2/3 


(3.24) 


Sarma et al. have not explicitly taken into account the effect of background /stray 
light intensity by introducing a parameter in the intensity equation. Instead, they 
proposed that one has to capture an image of unstressed model in the polariscope 
when polarizer and analyzer are kept in crossed position. The intensity values of 
this image are to be subtracted from all other images before any processing is 
done. Use of atan() function is recommended in this case. Figure 3.29 shows the set 
of experimentally recorded images. 



29 sequence of Four images recorded 

in Table 3 .7 for the ring under diametral compression ( loaa 3 

nm,F„= 11.2329 N/mm/fringe) 



3.6.6 Algorithm by Kihara 

Arrangements suggested by Kihara [1] are summarised in Tabie 3.8. 


uu 


Table 3.8. Optical arrangements for Kihara. 
Kihara 
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m 
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II 

CO 
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kI2 

nIA 

I 
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0 
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0 
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2^ .7^ 
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L 2 2 
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0 

mi 
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2 


-^/4 

m\ 

■ ; -nlA 


iS • oS . 9 

cos'^— + sm — sm 20 
[22. 


m\ 

mA 


» oS 7 

sm — cos‘^20 
2 



For the optical arrangement used by Kihara, the expressions for dete rminin g 
isoclinic and fractional retardation are 



1 

2 


tan’ 


V.^3 



1 _ifsin^sin2^'l 

-tan — 

2 l^sin<?cos20j 


for sin5 ^ 0 


(3.25) 


Use of atan() function is recommended to calculate theta from the above equation. 
Though the method proposed by Kihara is quite different, ultimately die method reduces 
to the one given by Patterson and Wang in view of the similarity in the intensity 
equations used for theta calculation. Figure 3.30 shows the set of experimentally 
recortted im^es. 








for the method by Kihara, as given 
I dia. 80 mm, thickness 5.05 mm, F< 


Fig 3.30 Sequence of first six images recorded experimentally 
Table 3.8 for the ring under diametral compression ( load 503 I 
1 1 .2329 N/mm/fringe ) • 


Fig 3.30 Sequence of last two images recorded experimentally for the method by Kihara, as ^yen in 
Table 3.8 for the ring under diametral compression ( load 503 N, dia. 80 mm, thickness 5.05 mm, = 

1 1 .2329 N/mm/fringe ) 

3.6.7 Comparison of Results obtained by Circular and Plane Polariscope 
Algorithms 

Figures 3.31 & 3.32 show the isoclinics obtained by using theoretically generated 
images. Figure 3.31a & b shows the isoclinic plots obtained for Patterson & Wang and 
Ajovalasit et al. methods. Figure 3.31c & d shows the isoclinic plots obtained for 
Ramesh «feMangal and Brown & Sullivan, methods. Figure 3.3 le & f shows the 
isoclinic plots obtained for Chen &Lin and Dupre et al. methods. Figure 3.32a & b 
shows the isoclinic plots obtained for Sarma et al. and Kihara. Methods. It is clear that 
all the algorithms give good results theoretically. But experimentally this is not the case. 
Figures 3.33 & 3.34 show the isoclinics obtained by using the experimentally recorded 
images. Figure 3.33a & b shows the isoclinic plots obtained for Patterson & Wang and 
Ajovalasit et al. methods. Figure 3.33c & d shows die isoclinic plots obtained for 
Ramesh &Mangal and Brown & Sullivan, methods. Figure 3.33e & f shows the 
isoclinic plots obtained for Chen &Lin without adjusiment and Chen &Lin wifli 
adjustment. Figure 3.34a & b shows the isoclinic plots obtained for Dupre et al. and 
Sarma et al. methods and Fig. 3.34e shows the isoclinic plots obtained for Kihara’s. 

method. 


Fig 3.31 Isoclinic plots obtained for the theoretically generated images of ring under diametral 
compression, a by Patterson Sc Wang b by Ajovalasit et al. c by Ramesh & Manga! d by Brown 
& Sullivan e by Chen & Lin f by Dupre et al. 






Fig 3.32 Isoclinic plots obtained for the theoretically generated images of ring under diametral 
compression, a by Sarma et al. b by Kihara 








Fig3.33.sooWop,o.ob«™dfo,^e^n«y^^^ 



Fig 3.34 Isoclinic plots obtained for the experimentally recorded images of ring under diametral 
compression, a by Dupre et al. b by Sarma et aL c by Kihara 

Figures 3.35 and 3.36 show graphs for the experimentally obtained data at lines y/R = 
0.5 and y/R = 0.25 respectively. The experimentally calculated theta is compared with 
the theta calculated theoretically 
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Fig. 3.36.C Comparison of Experimenatal theta with theoretical theta for line aty/R =0.25 for Chen & Lin, 

Sarma et. al. & theory 

Figure 3.35 a shows the graph obtained from the methods of Ramesh & Mangai, 
Patterson & Wang, Ajovalasit et. al. and dieory for the line at y/R = 0.5 It shows that the 
results obtained by the methodologies of Patterson & Wang and Ajovalasit are having 
sharp hikes. Thus the results obtained are not good. The results obtained by Ajovalasit 
et. al. are having lesser hikes than Patterson and Wang. This is because of error 
compensation by using left and right circularly polarized light beir^ used for model 
illumination. But it is clear that error is not fully compensated. Results by Ramesh and 
Mangai are closer to the theoretical line. Figure 3.35b shows the graphs obtained from 
the methods of Brown & Sullivan, Kihara, Dupre et. al. and theory for the line at y/R = 
0.5. It shows that the results obtained by Kihara’s technique are having sharp hikes, the 
data follows the trend of the theoretical line in between the hikes. However the results 
obtained are far from satisfactory. The results obtained by Dupre et. al. ^e not correct. 
Results by Brown and Sullivan are closer to the theoretical line. Figure 3.35 c shows the 
graphs obtained from the methods of Chen & Lin (without adjustment and with 
adjustment), Sarma et. al. and theory for the line at y/R = 0.5. It shows that the results 




obtained by Sarma etal. technique are having sharp hikes. Thus the results obtained are 
not good at some portions and on some portions they are close to the theoretical data.. 
The results obtained by Chen & Lin with adjustment are good and follow the theoretical 
curve closely. Chen and Lin without any adjustment is not good at some portions and 
are good in some other portions. Along line at y/R = 0.5 theta is not beyond the range — 
7 d\ to +7 z/ 4 and hence the utility of these algorithms for a generic situation is not fully 
varified. 

Hence the study was done for the line at y/R = 0.25. Along this line, the theoretical 
value of ^lies in the range of —nil to +7d2. Here the results are s imil ar to what has been 
observed previously. However none of the algorithms give correct value of ^beyond — 
tiJA to range. 

3.7. Closure 

Study of intensity patterns is carried out for the understanding of numerical inaccuracies 
involved in the isoclinic parameter calculations. The theoretical possibility of arriving at 
the correct and unwrapped value of isochnic parameter by using isochromatic parameter 
is considered. The quality of isoclinic data in circular polariscope for different loads is 
studied. Using isoclinics from low loads in circular polariscope is recommended. It has 
been reconfirmed that plane polariscope give a comparatively better isoclinic parameter 
even at higher loads as compared to the circular polariscope. 

From the study conducted in this chapter it is reaffirmed that the results obtained 
in plane polariscope are better than in circular polariscope. In circular polariscope, the 
method by Ajovalasit et. al. is found to be better as compared to the method by 
Patterson & Wang. Even though theoretically the principle stress directions remain 
constant with respect to the load, experimentally they are distorted at higher loads. In 
plane polariscope techniques, Ramesh & Mangal, Brown & Sullivan give better results. 
Method of Chen & Lin also gives good results but an adjustment is required. Methods 
of Dupre et. al. & Sarma et. al. are not giving satisfactory results. Besides this the theta 
obtained in all algorithms is wrapped and is in the range of -nJA to +nlA, fiius 

Unwrapping is still required from these algorithms. 


CHAPTER 4 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE 

WORK 


4.1 Conclusion 

The main aim of this thesis is to understand the nature and evaluation of isoclinic 
parameter by PST mediods. Methods to improve theta by quadrant method helped in 
correcting the phase map only for particular cases. The use of logical NOT has showed 
a way for correcting phase map for any generic problem. However, the mediod demands 
identifying the zones for correction interactively. 

It is established in the thesis that the orientation of model slow and fast axis has 
a significant effect on the final result of phase map. The study has showed that one can 
get a complementary phase map either by changing the direction of loading or by just 
swapping the intensities h & I 5 and U 8c h in the calculations. This has led to another 
study wherein the imwrapping of 9 based on intensity information has been carried out. 
However, due to interaction of <5 as well as ^ , it was not possible to identify whether 
one has model fast axis or slow axis at the point of interest based on intensity 
information ^one. A method to get 6 correctly in the range -Jiil to + -iri2 has been 
developed provided 5 is known for the problem a priori. The method has been shown to 
work based on theoretical data. Further refinement is required to apply it for practical 
situations. 

A comparative study on the evaluation of ^by various phase shifting methods in 
circular polariscope, plane polariscope and polarization stepping method have been 
carried out. The study showed that, in circular polariscope, the method by Ajovalasit et. 
al. is found to be better as compared to the method by Patterson & Wang. Even though 
theoretically the principle stress directions remain constant with respect to the load, 
experimentally they are distorted at higher loads. In plane polariscope techniques. 
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Ramesh & Mangal, Brown & Sullivan give better results. Method of Chen & Lin also 
gives good results but an adjustment is required. Methods of Dupre et. al. & Sarma et. 
al. are not giving satisfactory results. Besides this the theta obtained in all algorithms is 
wrapped and is in the range of-jdA to +7dA, thus unwrapping is still required from these 
algorithms. 

4.2.Suggestion for Future Work 

1. Theoretically isoclinic parameter does not change with change in load. But 
experimentally in circular polariscope isoclinic get distorted at higher loads as 
compared to lower loads under similar conditions. It will be useful to understand the 
exact reason for this anomaly. Note that in plane polariscope, isoclinics are good 
even at higher loads. 

2. It has been found that background intensity and incident light plays an important 
role in the accuracy of results. Special efforts can be taken to minimi ze the 
background intensity and to have a uniform incident intensity over the model. 

3 . An interpolation — extrapolation technique can be developed to unwrapp theta based 
on monitoring the phase changes while processing 

4. Special attention should be given to develop the equations and methods which will 

give ^directly in the range oi-Tdl to +nl2. 
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